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The VT— symmetric differential Schrodinger equation Hip =- Eip with the operator 

O 

a > it is re-arranged as a linear algebraic diagonalization. With rigorous proof, our 






H = H(x) =p 2 + ax 4 + iPx 3 + cx 2 + i5x = H*(—x) on L 2 (— oo, oo) is studied. At 



non- variational construction of bound states offers an infinite-dimensional analogue to 
the recent finite-dimensional quasi-exact solution available at the less common a < 0. 
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1 Introduction 

The study of the general parity-breaking anharmonic oscillators 

V(x) = ax 4 + bx 3 + cx 2 + dx (1) 

has a colourful and inspiring history [|IJ. Its latest turn came with the recent letter by 
Bender and Boettcher || who discovered that after a partial restoration of symmetry 
in the complex plane, 

V(x) = V*(-x), (2) 

these potentials may become solvable quasi-exactly ||. The extensive numerical ex- 
periments indicate that all the spectrum of energies E in the similar potentials is real, 
bounded and discrete. According to the conjecture by Daniel Bessis || this puzzling 
observation might be a straightforward mathematical consequence of the "weakened 
hermiticity" (Q) . With motivations ranging from field theory || and nuclear structure 
|| up to solid state physics |7j this hypothesis finds its further support in a few explicit 
analytic || and numerical |J constructions and semi-classical |10| or perturbative [ll[ 



arguments as well as in several available rigorous mathematical proofs pj 

In the particular model (|1|) the condition (|2|) means that the couplings a and c 
remain real while their partners b = i f3 and d = iS are purely imaginary. The atten- 
tion of paper [[J was solely paid to the negative values of the asymptotically dominant 
coupling a because of the related finite-dimensional reducibility and subsequent par- 
tial solvability of the bound state problem at certain special couplings and energies. 
In the present paper we intend to complement and complete the latter study by a 
parallel analysis of its "less solvable" alternative with a > 0. One has to imagine 
that in spite of the manifest non-hermiticity of the related Hamiltonian the procedure 



of quantization may be kept equally well defined at any sign of a [13 j. Extensive 



discussions of this point date back to the famous Dyson's argument (|I3|; cf. the very 



recent summary in [|15]). The same or similar formalism covers even the limiting case 



with the vanishing a = at any complex b ^ p6 . 

The main reason for the conventional a > is the simplicity of its physical inter- 
pretation. In an unquantized world the minus sign of a would mean that the particle 



can disappear and return from infinity in a finite time. In comparison, the asymptot- 
ically real and growing V(x) admits a more immediate intuitive understanding. Its 
choice weakens the impact of the unusual invariance (0) which, in certain applica- 



tions, mimicks the combined effect of parity and time reversal VT ||10| . With a > 
we also get in a closer contact with the already existing calculations |T7J and with the 
re-summations of perturbation series, say, in terms of the so called Hill determinants 



I8| , |T9f or analytic pD| and matrix |2T| continued fractions. The new, non-Hermitian 
options (3 7^ and 5^0 open the new perspectives. 

We intend to show that the positivity of a which excludes the quasi-exact solv- 
ability need not contradict an efficient linear algebraic description of bound states. 
We shall re-write the differential Schrodinger equation in an equivalent matrix form. 
Although its dimension remains infinite, its structure and derivation will parallel its 
Hermitian Hill-determinant predecessors characterized globally by a loss of their her- 
miticity (pars pro toto, the reader may consult the review |£J). For our present 
non-Hermitian interaction (P + (0) such a loss is much less harmful. 

2 Non-terminating recurrences at a > 

Forces ([I]) with the "weak" symmetry (|2[) enter the differential Schrodinger equation 

( d 2 \ 

— -— + ax 4 + i (3 x 3 + ex 2 + iS x J ip(x) = E ip(x), x £ (—00, 00) (3) 

which has the two independent asymptotic solutions 



^j {±) {x) 



exp 



3 2 

X X 



u— + v— + 0(x) 



u = ±y/a=£0, v = i—. (4) 

2m 



This explains the difference between a < and a > 0. In the former case we may 
move the real axis downwards in the complex plane, x — r — irj, rj > 0, r £ (—00, 00). 
Whenever we pick up 77 > — /5/4w |a|, we discover that the asymptotic solution ip(~'(x) 
remains integrable at both the ends of the real axis r. Some of the (necessarily, 
analytic) bound states ^jj( exact ) ( x ) may (and do) acquire an elementary form of an 
exponential-times-polynomial product for a < ||. 



Let us now consider a positive coupling a > and re-scale its value to a = 1 
for simplicity. The general solution of our Schrodinger equation (§) will have the 
form ip( gen \x) = n±ip <y+ \x) + v±ip(~'($). It may only be made compatible with the 
required asymptotic decrease near x —* ±00 by the sign-of-x-dependent choice of its 

parameters, 

^(pfc»«)( x ) = I / + ^(-)(a;), x>l, 

( 5 ) 
tfj(p h y s )( x ) = /i_^(+)(x), x<-l. 

In contrast to the preceding case, the exact solution cannot be constructed as a 

product of an exponential with a polynomial anymore. One has to resort to the next 

eligible possibility, say, 

00 
^ ansatz \x) = e~ sx2 Y, h n {ix) n . (6) 

n=0 

This is a manifestly VT— invariant infinite-series ansatz. Abbreviating ix = y we 
derive the recurrences 

A n h n+2 + C n h n + 5 h n -\ + h n _ 2 - (3 h n - 3 + h n _ A = 0. (7) 

All its coefficients are real, A n = (n + l)(n + 2), C n = 4sn + 2s — E and 9 = 4s 2 — c. 
By construction, the set (H) is equivalent to our differential equation (|3|). At all 
the tentative energies E it defines the coefficients h n from an input pair ho and hi. 
We have to determine these parameters via a fit of (§) to the appropriate boundary 
conditions 

^ anmtz \X R )=0 = ^ ansatz \-X L ), X B >1, X L >1. (8) 

In comparison with the other numerical methods of solution of the bound state prob- 
lems with symmetry (0) || such a recurrently specified recipe does not look any 
superior, especially because it requires a cumbersome numerical limiting transition 
Xr,l — > 00. A deeper insight and simplifications are asked for. 

3 The asymptotics of coefficients /i„ at n > 1 

Recurrences (|7|) form a linear difference equation of the sixth order. One may recall 
the standard theory of its solution [53] as well as its immediate application to quartic 



oscillators p4|, p5| . At the large indices n, the sextuplet of the independent asymptotic 



solutions h n acquires the general Birkhoff form as presented, say, in ref. [26]. All 
these solutions decrease as h n ~ 0(n~ n ^ 3 ) at least. For our present purposes, they 
may easily be re-derived as follows. 

Firstly, in the leading-order approximation, we replace equation ([?]) by the mere 
two-term dominant relation between h n+ 2 and /i„_4. This inspires us to change vari- 
ables h n — > g n and we re-write all the six independent solutions in the same compact 
form 

^) = ( 3 i/3)nr(l + n/3) ' P=l,2,.-,6. (9) 

The p— dependent complex parameter X(p) = exp[z(2p — l)7r/6] characterizes the dom- 
inant n— dependence of the separate solutions while the new functions or coefficients 
9n = 9n{p) vary more slowly with n. 

This confirms the linear independence of our six solutions but leaves their absolute 
values indistinguishable. In order to remove this degeneracy in size we re-introduce 
equation (|7|) in its amended, second-order asymptotic form 

9n + 2 ~ 9n-A = -^ 9n ~ ^ 9n-3 + O [^ J . (10) 

With most of the components of the Stirling formula still hidden within the error 
term, the smallness of the ratio 1/n 1 ^ 3 enables us to infer that 

g n = e 7n^ +0 („i/3 )} n>1 (n) 

The complex exponent 7 = 7(2?) = s\ 4 (p) — /3 \(p)/4 depends on p. An elementary 
trigonometry gives the explicit formulae 

Re 7 (1) = Re 7(6) = ~P - |, Re 7 (2) = Re 7 (5) = s, 

Re 7 (3) = Re 7 (4) = ^-|. (12) 

In combination with eq. @ this already implies that the radius of convergence of our 
Taylor series (|6|) is infinite. The function ^p( ansatz ) ( x ) is unique and well defined at any 
complex x. Its shape is fully determined by the energy E and by a not yet specified 
choice of the two initial complex coefficients ho and h\. 



The second important consequence of identities (|T2| ) is that whenever we satisfy 
the condition 

-4 (13) 

the general solution h n = J2p=iG p h n (p) itself will be asymptotically dominated by 
its two most quickly growing components, 

h n = G 2 h n {2) + G s h n {5), n »1. (14) 

In this sense we are free to set G\ — G% — G4 — G$ — in the asymptotic domain 
of n ^> 1. Each choice of the energy E and initial h and hi will only generate a 
different, x— and n— independent pair of coefficients G2 and G 5 . 

4 The asymptotics of <ip( ansatz ) (V) at \x\ ^ 1 



Equation (14) is a key to our forthcoming replacement of the numerically awkward 
boundary conditions (§) by the much more natural approximative truncation of recur- 
rences (|7|). We shall parallel the Hermitian construction of ref. |22| and try to bracket 
the exact energy between its upper and lower estimates with E ^ E(physical). Un- 
der such an assumption our infinite series ■ l p( ansatz ) ( x ) as defined by equation (J^j will 
always exhibit an exponential asymptotic growth as described quantitatively by eq. 
(|j) above. This means that we shall exempt the possible lucky guess of the exact 
energy in its full precision as never relevant in any step of our forthcoming consider- 
ations. Such a very formal point of view does not contradict the underlying physical 
intuition since boundary conditions (§) are approximative. One has to move to the 
limit Xr,l — > oo in principle. 

The most important immediate consequence of our "bracketing" interpretation of 
boundary conditions is that at the large absolute values of the coordinate \x\ 3> 1 the 
first N exponentially small components 0(e~ sx ) may safely be ignored as irrelevant. 
We may also insert (|j) and (|T4"D in ^p( ansatz ) ( x ) ~ exp(— sx 2 ) J2™=n+i h n (ix) n with 



iV > 1 and get 

b (ansatz) ( )r ^ -sx 2 y> G 2 \ n (2) g n (2) + G 5 \ n (5) fl n (5) 



The validity of this formula is a strict consequence of the specific constraint ([13] 
imposed (say, from now on) upon the admissible quasi-variational parameter s. 
Once we split ip( ansatz ) ( x ) = ^p( ansatz ^G 2 , G$,x) in its two components 



v> 



(ansatz) 



(G 2 ,0,x)~G 2 e- sx J2 

n=N+l 



K —x) exp 



7 (2) n 2 '* + 0{n x ' z ) 



(3V3)nr(i + n /3) 



ip ( - ansatz) {0,G 5 ,x)~G 5 , 



E 

n=N+l 



x exp 



7 (5)n 2 / 3 + C(n 1 / 3 ) 



(3V3)»r(i + n/3) 



we may apply the rule e z ~ (1 + z/t)*, t ^> 1 in the error term and get 



V> ( 



ansatz) 



(GiA-y) 



exp[—sy z 



and 



^ mss(2 )(0,G 5l 2/) 



G2 J + i(3 1/3 ) n r(i + »/3) 



G5 „=? + i(3 1/3 )"r(l + n/3) 



!/• 



1 + 



ATi/3 



V- 



1 + 



ATi/s 



exp(— sy) 

This is valid at all the large arguments y. Along the positive semi-axis y ^> 1, 
both the right-hand-side summands are real and positive. They sum up to the same 
function exp[y 3 /3 + 0(y 2 )}. This is a consequence of the approximation of the sum 
by an integral and its subsequent evaluation by means of the saddle-point method. 
The same trick was used by Hautot, in similar context, for the V— symmetric and 



Hermitian anharmonic oscillators [27]. 

In contrast to the Hautot's resulting one-term estimates of ip, the present asym- 
metric, VT— invariant construction leads to the more general two-term asymptotic 
estimate 

^ ansatz) (G 2 , G 5 , x) ~ G 2 exp[-x 3 /3 + 0{x 2 )} + G 5 exp[a; 3 /3 + 0(x% \x\ > 1. 

As long as we deal with the holomorphic function of x, this estimate may be analyti- 
cally continued off the real axis of x. Near both the ends of the real line and within 
the asymptotic wedges |Im x|/|Re x\ < tan7r/6 we simply have the rules 

^ ansatz \G 2 ,G 5 ,x) ~ G 2 exp[-x 3 /3 + 0(x 2 )}, Re x < -X L < -1 (15) 

and 

^ ansatz \G 2 ,G 5 ,x)~G 5 exp[x 3 /3 + 0(x% 1< X R < Re x. (16) 

They are fully compatible with formula (151) since a — 1. 



5 The matrix form of the Hamiltonian 

Our complex differential Schrodinger equation (|3|) becomes asymptotically real, in the 
leading-order approximation at least. In a suitable normalization the wave functions 
^p(ansatz) ^ ma y ^ e mac [ e asymptotically real as well. Near infinity they will obey 
the standard Sturm Liouville oscillation theorems f28|. In particular, after a small 
decrease of the tentative energy parameter E > E(physical) the asymptotic nodal 
zero Xr or — X^ originating in one of our boundary conditions (^) will move towards 
infinity p|. 

This may be re-phrased as follows. At a more or less correct physical real pair 
h Q = p cosC and h\= p sinC with ( G (0, 2tt) and with the convenient normalization 
p = 1 a small change of the energy E somewhere near its correct physical value 
Eq « E(physical) will cause a sudden change of the sign of the asymptotically growing 
exponentials ( |i5|) and fll6| ) at some Co ~ ({physical). This may be re-read as a doublet 
of conditions 

G 2 = G 2 (E , Co) = 0, G 5 = G 5 (E , Co) = 0. (17) 

In the limit iV — > oo of vanishing corrections, these two requirements may be re- 
interpreted as a rigorous re-incarnation of our original physical asymptotic boundary 
conditions (^). The conclusion has several important consequences. Firstly, at a fixed 
JV>1 we may define 

X N (p) exp[ 7 (p)iV 2 /3] 
Jp-^P ( 3 i/3)Jv r ( 1 + JV / 3 ) > P-4&- 

Functions / p = / P (-E, Co) differ from their sign-changing predecessors G p = G P (E, Co) 
just by a constant factor near E , f p (E, Co) ~ F p ■ (E — E ). We may write 

/^ « (F 2 + F 5 )(E - E ) + 0[(E - E )% 
(N + 3y/ 3 h N+1 « [F 2 A(2) + F 5 A(5)](£ - £ ) + 0[(E - E Q f] 

due to equation (Q). This formula connects the two functions G 2 , G5 with the two 
neighboring Taylor coefficients h N = h^(E Xo) an d h N+1 = h N+1 (E ,(o) near the 
physical E and Co by an easily invertible regular mapping. This means that the 



implicit algebraic boundary conditions ( |T7| ) are strictly equivalent to the fully explicit 
requirements 



M^o,Co) = o, 



h N+1 (E , Co) = 0, iV>l. 



By construction, this becomes an exact physical bound-state condition in the limit 
N — > oo. At the finite iV ^> 1 its appeal lies in its change-of-sign character. This 
need not make equation fll8l) immediately suitable for computations but once we fix 
N = N ^ 1, E = E , ( = ( and insert the zeros ( |TE| ) in our recurrences (|7|), we 
arrive at the truncated square-matrix equation 



/ 



Co A 

S C X Ay 

e 5 ••. ■•. 

-p e ••. 

1 -{3 ••• 



V 



1 



-(3 



'•• ^iV-3 



5 C N -i J 



\ 


( h \ 




h 




h 2 




h-N-3 




hN-2 


) 


y h,N~i j 



0. 



(19) 



This is our main result. As long as C n = 4sn + 2s — E, the energy enters just the 
main diagonal and we may determine all its approximate low-lying values Eq by the 
routine iV x iV— dimensional diagonalization. 

6 Discussion 



6.1 Illustrative numerical tests 

The smallest matrix in equation ([19]) which contains all the couplings has dimension 
N = 5. It is quite surprising that such a drastic simplification leads to the mere 5% or 
6% error in the ground-state energy. Together with the equally pleasant quick increase 
of precision with the growing N, this is illustrated in Table 1. Table 2 shows where 
the numerical application of the present approach can find its natural limitations. We 
observe a steady decrease of precision at the higher excitations. 



With s = 2 and a 



(3 = 5 = 1, both Tables were computed in MAPLE 



[ 30f] . This language keeps the possible loss of precision under a careful control. This 
implies the growth of the computing time at the higher dimensions. Still, the very 
quick actual numerical rate of convergence enabled us to compute all our examples 
on a current PC in a couple of minutes. 



6.2 Determinantal formulae for the Taylor coefficients 

With real ho and hi our wave functions ip( ansatz ) (x) are composed of the spacially sym- 
metric real part and spacially antisymmetric imaginary part. Such a normalization 
fixes the phases of the complex constants G 2 and G5 accordingly, i.e., via equation 
(PI)- This clarifies the structure of the asymptotics of the wave functions. 

Polynomial approximants of the Taylor series (Q offer a reliable picture of i]){x) 
in a broad vicinity of the origin. We may recall recurrences (^) and reveal that the 
ho— and h\— dependence of any coefficient h n is linear, 



h n = h a n + hiio r , 



CTq = OOl = 1, o~i = LUq = 0. 



All three sequences h n , a n and u n satisfy the same recurrences. As long as <j\ = 
and a>o = we may omit the second or first column from equation (|l^) in the latter 
two respective cases. In terms of the (m + 1)— dimensional matrices 



/ 



Co 


Ao 












s 




c 2 


-4i 



A 2 








1 




c 3 













-/? 









A m -2 






1 




8 


C m _i 

s 



r 


A m -\ 





\ 



and 



a 



V 






^4o 












Ci 





^i 










<5 


<?2 

<5 



C 3 


^2 








1 




5 







A- m -2 








-P 




^m— 1 





A m -1 






1 




s 


L"m 






/ 



we may re-write not only the recurrences themselves but also their unexpectedly 
compact solution 



C"n+1 



:-i: 



detS 



n-l 



detfi 



k>n+l 



n-l 



n 



1,2, 



(20) 



nl(n+l)V ""^ v ' n\(n + l)\ 

We need to know just the correct physical values of the three variable parameters 
(viz., the norm p = Jh^ + h 2 , the ratio hi/ha = tan£ and the physical energy E) 
in order to be able to define our physical wave function ip(x) completely in terms of 
these closed formulae. 



6.3 Alternative ansatzs and constructions 

We have shown that the Taylor-series ansatz (§) mediates a useful transition from 
differential equation (|3]) to the difference equation (0), followed by its further replace- 
ment by our final matrix Schrodinger equation ([TJ^). In this context it is important 
to mention that our choice of the initial form of ansatz (^ is by far not unique. 

A nice example of an alternative expansion may be found in ref. |31| where the 
Hill-determinant study of the symmetric potentials V(x) = x 2 + A u(x) with the non- 
polynomial anharmonicity u(x) = x 2 /(l + g x 2 ) via the series of the form (|]) has been 
rendered possible by the use of the Taylor series in powers of the "adapted" variable 
u(x). Sophisticated versions of the latter trick move the (complex) singularities off 
the physical domain of convergence and their active use in physics dates back to Jaffe 
32|j at least. They may even help us to deal with relativistic corrections |33| etc. 



10 



Unfortunately, an application of the changes of variables to asymmetric potentials 
is not without its specific difficulties. Efficient methods of their suppression have been 
suggested, therefore, in our older paper |34j and, recently, by Bay et al [p5fl . Most 



often, one employs the two independent separate ansatzs (one for each half-axis) and 
matches the wave functions, say, in the origin. 



In the latter comparison, the method of paper |25| is most straightforward. It is 
based simply on an introduction of the second free parameter (G or V in the original 
notation). Even from the very numerical point of view, the essence of the algorithm 
of Bay et al remains purely iterative, therefore. 



The more algebraic method of ref. [34] works directly with the matched, "doubly 



infinite" sparse matrices. Although the algorithm itself is already fully algebraized, its 
universality seems redundant for our present purposes. Indeed, the VT— symmetric 
forces (|]) are only composed of the real part which is spatially symmetric and of 
the non-vanishing imaginary part which is spatially antisymmetric. This additional 
information is well reflected and used by our present non-matching approach. 

We may summarize that we were able to preserve a maximal similarity of our 
"new Hill determinants" to their current Hermitian predecessors (cf., e.g., ref. |nj). 
Moreover, in a way completing the parallel studies of the other VT symmetric po- 
tentials, the very specific form of their spatial asymmetry proved again "extremely 
weak" from the purely methodical point of view. We re-confirmed that its simplifying 
role strongly resembles the role of the usual V— symmetry (i.e., parity), so useful in 
many parts of the current textbook quantum mechanics. 
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Table 1. 

The N— dependence of energies. 



N 


E 


E 1 


5 


1.793 


7.547 


6 


1.823 


5.868 


7 


1.634 


5.856 


8 


1.673 


5.138 


9 


1.627 


5.162 


10 


1.658 


4.922 


15 


1.693 


5.106 


20 


1.692 


5.126 


21 


1.691 


5.124 


22 


1.692 


5.123 


23 


1.692 


5.123 


24 


1.692 


5.123 


25 


1.692 


5.123 
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Table 2. The growth of precision with dimension N. 











energies 


E n 








n 





1 


2 


3 


4 


5 


6 


7 


N 


















15 


1.69347 


5.106 


9.152 


13.043 


17.817 


23.89 


31.26 


41.55 


20 


1.691638 


5.12559 


9.2800 


14.050 


19.244 


— 


32.35 


— 


25 


1.691579 


5.123441 


9.25812 


13.8689 


18.7925 


24.265 


30.039 


37.97 


30 


1.691590 


5.123614 


9.26174 


13.8826 


18.8922 


24.262 


29.726 


34.67 


35 


1.691590 


5.123579 


9.26151 


13.8793 


18.8838 


24.220 


29.860 


35.85 
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